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Q i ABSTRACT 

o 

; We describe the results of numerical simulations of the dynamics of the bound- 

ary layer (BL) between the accretion disk and the surface of a non-magnetic 
white dwarf (WD) for different viscosities which corresponds to different stages 
for dwarf novae burst cycles. The simulations cover the inner part of the accre- 
(N ■ tion disk, the BL, and the upper atmosphere of the star. The high viscosity case, 

which corresponds to a dwarf nova in outburst, shows a optically thick BL which 
after one Keplerian rotation period (t^=19s) extends more than 30 degrees to 
either side of the disk plane. The BL is optically thick and thus occludes part of 
^ ■ the star. The low viscosity case, which corresponds to a dwarf nova in quiescense, 

also shows a BL, but it is optically thin. 
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1. Introduction 

A cataclysmic variable is a close binary system where gas flows from a Roche-lobe 
overflowing low- mass star via an accretion disk to the surface of a WD (Warner 1995). 
This flow releases gravitational energy at a rate of L acc = GM*M / R*. In the disk, the gas 
describes almost circular Keplerian orbits with vk < y 'GM^j 1 R*. If viscosity is present in 
the disk, viscous shear forces between the differentially rotating disk annuli cause friction 
which dissipates energy and turns the orbital energy into heat. The heat is radiated away 
and the resulting energy loss causes the gas to fall into a lower orbit (Shakura & Sunyaev 
1973). This is the mechanism behind the accretion process. 

In many accreting WD systems, such as dwarf novae and novae-like variables, the mag- 
netic field is weak and the disk extends all the way to the WD surface. In order to be 
accreted, the matter in the accretion disk must lose all of its remaining rotational kinetic 
energy and decelerate from a nearly Keplerian velocity to the much slower rotational velocity 
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of the surface. This happens at a rate of ~ \Mv\ \GM*M / 1 which means that about 
one half of the gravitational energy is released in the BL between the disk and the WD as the 
rotating matter settles down on the surface while the other half of the gravitational energy 
drives the disk luminosity (Lynden-Bell & Pringle 1974). Here the detailed transition of the 
gas flow determines the relative fractions of how much of the kinetic energy is dissipated as 
heat and radiated away (Pringle 1981), radiated or advected into the surface of the WD, and 
how much is converted into winds, and how much goes into spinning up the WD and heating 
it up. Additional heat will be released from compressing the atmosphere. The details of 
radiated spectrum, therefore, depends on the structure and dynamics of the BL! 

The transition region has typically been treated by a BL model, where it has been 
assumed that the averaging of the vertical structure over a disk scale height, H (see Pringle 
(1981)), which is employed in accretion disk studies (Shakura & Sunyaev 1973) can be 
extrapolated to the surface. This reduces the description of the BL to a one dimensional 
problem (Pringle 1981; Meyer & Meyer-Hofmeister 1982; Popham & Narayan 1995; Collins 
et al. 1998a,b). Ferland et al. (1982) pointed out that this radial extrapolation might be 
invalid. Inogamov & Sunyeav (1999) (IS99) presented the first analytic treatment of the 
latitudinal aspect of the BL. Averaging over one pressure scale height, h = Pdr/dP, in the 
radial direction Piro & Bildsten (2004) (PB04) calculated the vertical structure of the BL 
for a WD and found that the BL could spread towards the poles of the WD. Early numerical 
calculations used an axisymetric two-dimensional configuration (Robertson & Frank 1986; 
Kley 1991), but since their simulations predated the papers of IS99 and PB04, they did not 
use enough resolution in the radial direction to capture the dynamics of the BL. 

The objective of the paper is to determine the structural and dynamical dependence 
of the BL as parametrized by the viscosity. In section 2, we describe our well- resolved 
simulations that accurately capture the dynamics of the BL. In section 3 we discuss the 
multidimensional flows that arise in the BL and present the differences between flows with 
different accretion rates. 

2. Numerical setup and model 

Any numerical simulation must resolve the characteristic scales of the object that is 
being simulated. The scale of the BL is equal to the star's atmospheric scale height in the 
radial direction and the disk scale height in the vertical direction. Furthermore a good model 
of the BL must include both the inner part of the disk and the outer layer of the atmosphere, 
so we include five atmospheric scale heights of the stellar atmosphere and three disk scale 
heights of the disk in our simulation. 
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To do this we use a spherical mesh with 384 ratioed zones in the radial direction and 
128 ratioed zones in the azimuthal range spanning to 30 degrees from the disk plane. This 
is sufficient to resolve a BL of a few atmospheric scale heights with more than 100 zones for 
a very hot star with a temperature of 3.0 x 10 5 K. The disk has a temperature at the base 
of To = 1.0 x 10 6 K and it is surrounded by a halo, which is 50 times hotter than To. The 
star, disk and halo are tagged by two different species tracers in the simulation. The toroidal 
component is assumed to be symmetric around the axis and a reflective boundary condition is 
imposed at the midplane. This approximation is justified because the gravitational distortion 
of the secondary is negligible this close to the primary. The equilibrium disk and halo model 
setup is described in Balsara (2004) which builds on models described in Suchkov et al. 
(1994) and Matsumoto et al. (1996). 

We solve the compressible Navier-Stokes equations as given by Mihalas & Mihalas (1984) 
using the RIEMANN code. The inclusion of viscous terms in the momentum equation also re- 
quires a self-consistent introduction of viscous dissipation into the energy equation. We follow 
the prescription of Mihalas & Mihalas (1984) in our treatment of viscous energy dissipation. 
The spatially and temporally second order agorithms in RIEMANN have been described in Roe 
& Balsara (1996); Balsara (1998a,b); Balsara & Spicer (1999a,b) and Balsara (2004) and use 
many ideas from higher order WENO schemes (see Jiang & Shu (1996) and Balsara & Shu 
(2000)) to reduce dissipation. 

The matter in the model is subject to the central gravitational field of the underlying 
WD, which has a radius and mass of = 9 x 10 s cm and M* = O.6M respectively (same 
as PB04). The viscosity, which can be molecular, turbulent, magneto-rotational, radiative 
or a combination thereof, is described by a single parameter, a, as in Shakura & Sunyaev 
(1973), where v = ac s H (c s is the sound speed and H is the vertical scale height of the disk). 
In this letter we compare a = 0.1, which is representative of the disk in the outburst stage 
(Smak 1984; Mineshige & Oasaka 1985) with a = 0.001 which corresponds to a disk in the 
quiescent stage. The model uses an ideal gas (7 = 5/3) of a fully ionized solar composition 
(/i = 0.62g/mole) and assumes no radiation transport. 

3. Results and discussion 

Figs. 1 and 2 show the logarithm of the density with overlaid poloidal velocity vectors 
for the a — 0.1 and a = 0.001 cases. They show that a BL forms at all accretion rates thus 
corroborating the scenario of IS99 and PB04. The viscosity has a significant impact on the 
structure of the BL. 
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For a = 0.1, which corresponds to a dwarf novae in the outburst state, the infalling 
matter depresses the surface and forms a band of disk material that spreads out from the 
footpoint of the disk. Here matter keeps piling up as viscous shear forces prevent it from 
quickly moving to the sides. The shear forces between the rotating band and the underlying 
star spin up the much denser star matter while slowing down the band. As the band keeps 
evolving, friction dissipates energy which heats up the band. The pressure in the band 
increases due to both the increase of internal energy from viscous dissipation and the increase 
of density from accretion. At the disk-star interface, the accreting disk matter slows down 
to very sub-Keplerian velocities while the stellar atmosphere spins up. At higher latitudes, 
where the BL is not in contact with the accretion disk, the outer layers of the BL rotate 
faster than the inner layers. With decreasing latitude, the rotation in these outer layers of 
the BL blend in smoothly with the rotation of the disk. At lower latitudes, where the BL is 
in contact with the accretion disk, the pressure and the centrifugal force in the BL prevent 
the disk matter from accreting radially onto the star. Thus the pressure and the centrifugal 
forces in the BL force the disk matter to flow to higher latitudes where the pressure is lower. 
This further enhances the process of forming a multi-dimensional BL. 

In situations with high accretion rates, say a = 0.1 in Fig. 1, the ram pressure of the 
infalling matter can be sufficiently high to depress the surface of the star. Just as blowing 
down on the surface of a cup of coffee results in surface waves on the coffee, the depression 
on the surface of the star causes surface waves to be set up on the star. However, since 
the disk matter piles up too quickly, it soon overflows the surface wave and continues above 
it towards the pole at supersonic velocities. This raises the possibility that the poloidal 
flow is also susceptible to the rapid development of gravitational and/or Kelvin Helmholtz 
instabilities. Such instabilities have been studied in isolation by Alexakis et al. (2004). As the 
overflowing material spills over the surface wave, it creates a substantial amount of friction 
and heating at the point where it re-establishes strong contact with the star's surface. A 
hotspot develops at that point of contact, resulting in a reduced density, as can be seen 
starting at (r,z) = (0.87,0.23) in Fig. 1. This effect is not seen at lower accretion rates, see 
Fig. 2. 

The inclusion of opacity and radiative transport would cause part of the viscously dis- 
sipated energy to be radiated away, resulting in a decreased pressure in the BL. Since the 
pressure contributes to the spreading, it would cause a decreased spreading of the BL. The 
logical extreme of fully efficient radiative transfer corresponds to a situation where all the 
dissipated energy is instantly radiated away. This is equivalent to the exclusion of the vis- 
cous dissipation terms in the energy equation. We have carried out simulations in this limit 
also, and we find that the development of the BL is, nevertheless, a robust conclusion. Such 
simulations will be presented in a subsequent paper. 
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Eventually the accretion process should settle into thermal equilibrium as energy losses 
would cause the accreted band to sink and become part of the WD's atmosphere as fast as 
new matter accretes from above. Since we do not include radiative transport we do not trace 
this thermal equilibrium in the present simulation. However, for a = 0.1 (see Fig. 1) the 
BL extends more than 30 degrees away from the disk plane and reaches a height of 650km 
(~ 0.07-R*) after close to a Keplerian rotation period. The vertical extent of the BL is much 
larger than the extent of the inner disc as shown in Fig. 1 Consequently substantial heating 
of the inner disk by radiation from the BL may change both the density distribution and 
flow pattern near the inner edge of the disk (Calvet et al. 1992). In this paper we do not 
focus on the details of spectral modeling. Since Thomson scattering will always be present 
in an ionized plasma, we use it to provide a lower bound on the opacity. Here we note that 
the spreading layer would be opaque (r = — npdr > 2/3) to Thomson scattering up to 
at least thirty degrees away from the disk plane on either side. For a = 0.001 (see Fig. 2) 
the BL expands up to a thickness of about 110 km (~ 0.01/?*) after one Keplerian rotation. 
Here, matter simply flows towards the pole in a physically and optically thin layer that is 
approximately 50 km deep (~ 0.006-R*). 

Observations of the UV spectrum of U Geminorum during quiescense (Kiplinger et al. 
1991; Long et al. 1993; Froning et al. 2001; Szkody et al. 2002) show that the WD has a 
surface temperature of ~ 30000-38000K. However, additionally there are UV flux levels as 
well as metal absorption lines corresponding to solar abundances which can only be explained 
by the additional presence of a small hot (~ 57000K) and slowly decaying continuum source. 
N V and He II lines, which can only be formed in a hotter (T > 60000K) and optically 
thin region, are also observed during quiescense. We suggest that these observations may 
be explained by the transition from the optically and physically thick BL of the outburst 
stage (Fig. 1) to the optically and physically thin BL corresponding to the quiescent stage 
(Fig. 2). 

In conclusion, we have verified that BLs are multi-dimensional as suggested by IS99 
and PB04 and analyzed the multi-dimensional hydrodynamics of the BL for high and low 
accretion rates. These BLs cover a significant fraction of the surface of the star. For high 
rates of accretion the BL can be optically thick enough to occlude parts of the stellar surface. 
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Fig. 1. — The plot shows the outermost part of the star (to the left), the BL and the 
innermost part of the accretion disk (coming in from the right) for the a = 0.1 case. The 
baseline of the accretion disk corresponds to the bottom of the plot. The plot only shows the 
most interesting part of the entire computational domain. The colors show the logarithm 
of the density in cgs and the axes are in units of 10 9 cm. Velocity vectors in the poloidal 
direction are overlaid. The result corresponds to almost one Keplerian rotation period. 
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Fig. 2. — This plot shows the logarithm of the density in cgs for the a = 0.001 case at 
one Keplerian rotation. Notice the absense of waves on the stellar surface and the uniform 
thickness of the BL. 



